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A SIMPLIFIED METHOD FOR CALCULATING AEROELASTIC 

EFFECTS ON THE ROLL OF AIRCRAFT 

By John M. Hedgepeth, Pa^ll G. Waner, Jr . , 
and Robert J. Kell 


SIM4ARI 


An approximate linearized lifting-surface theory is used in conjunc- 
tion with structural influence coefficients to foimiulate a method for 
analyzing the aeroelastic behavior in roll of an aircraft. Rolling effec- 
tiveness and aileron-reversal speed are computed by the use of a Galerkin- 
type procedure. Results obtained for two example configurations by using 
this method are compared with the resialts obtained by using the more 
refined method of NACA TN 3067 . The agreement is excellent. 


INTRODUCTION 


In the design of modern high-speed aircraft, it is generally recog- 
nized that aeroelastic effects must be accomted for accurately. One 
method which should be capable of yieldi n g reliable predictions of the 
aeroelastic effects on the roll of supersonic aircraft has been presented 
in reference 1. This method, which males use of structviral influence 
coefficients to determine the distortions and lifting -surface theory to 
determine the airloads, involves, however, a considerable amount of com- 
putational labor. For this reason, some means for simplifying the com- 
putations without introducing an objectionable amount of error was sought. 
The purpose of this paper is to describe the restilting simplified method 
and to evaluate its accuracy. 

In this paper, attention is confined to the rolling problem. The 
actual aircraft configuration is left general, the only restriction being 
that the effects of chordwise bending are assumed to be negligible. Both 
subsonic and supersonic speeds can be treated by the method, but partic- 
ular attention is paid to the supersonic regime in the examples. 
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SmBOLS 


A parameter defined by eq^uation (15a-) 

B 0 ,Bp,B 5 parameters defined by equations (l 5 t>) 


Ciq,Ci^,Ci^ rolling-moment derivatives defined by equations (l 5 c) 


GJ 

L(y) 

M(y) 

Mo 

Ph 


elementary torsional stiffness of wing 

structural-twist influence function due to unit concentrated 
load at y-axis 

structural-twist influence function due to unit concentrated 
torque 

aerodynamic load per imit span, positive upward 

aerodynamic pitching moment about y-axis per unit span, 
positive in positive twist direction 

free-stream Mach number 

static pressure at altltiide 


PO 

Q(y) 


V 

a 

b 

c(y) 

c 




standard static pressure at sea level 

aerodynamic pitching moment about elastic axis per unit 
span, positive in positive twist direction 

free-stream velocity 

ratio of fuselage radius to exposed semlspan of wing 

total wing span, 2{a.l + 1) 

wing chord 

mean geometric chord 

aileron chord 

section lift coefficient, L(y)/qc 



NACA ™ 3370 


3 


Cm 


section pitclalng -moment coefficient about y-axis. 


M(y) /ic2 


e(y) 

k 

I 

P 

pb/27 

x,y 


P 

0(y) 

0i(y) 

5 

7 


section pitching-moment coefficient about elastic axis^ 

Q(y) /q.c2 

distance measured forward from y-axis to elastic axis, 
expressed as fraction of local chord 

amplitude of twist mode shape 

exposed wing semispan 

ro3-ling angular velocity, positive in right-hand sense 

tangent of wing-tip helix angle 

dynamic pressure 

coordinate system (see fig. 1) 

cotangent of Mach angle, - 1 

angle of twist of wing (see fig. 1) 

twist mode shape \ised in Galerkin-type procedure 

rolling effectiveness, (pb/2V')j>^(pb/2V)p 

aileron deflection (see fig. l) 
ratio of specific heats 


Subscripts : 


F 

R 

oe 

P 

rev 


flexible wing 
rigid wing 

effective aerodynamic coefficients due to twist 
parameters or aerodynamic coefficients due to unit pb/27 
aileron reversal 

parameters or aerodynamic coefficients due to unit aileron 
deflection 


6 
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parameters or aerodynamic coefficients due to unit twist 
shape 


Matrix notation: 


C J 
L J 
I I 

N 


square matrix 
row matrix 
column matrix 
diagonal matrix 


MALYSIS 


The analysis proceeds adong the same lines as that in reference 1 ; 
that is, the structiiral deformations are expressed in terms of the air- 
loads, the airloads are obtained, and then the two are combined to form- 
ulate the aeroelastic problem. 


Structural Deformations 

Consider the configuration shown in figxire 1 . If the effects of 
chordwise bending are assumed to be negligible, the only distortion of 
interest in this problem is the twist 9 (y) which can be expressed in 
terms of the section lift L(y) and the section moment about the y-axls 
M(y) as follows: 


I nl 

GL(y^Ti) L(ii) M(ti) dq (1) 

'^0 

where Gp(y,T)) and QM(y,Tl) are Influence functions which define the 
wing twist at y due to a unit concentrated load at the y-axis and a 
unit concentrated torque, respectively, at the spanwise station t]. As 
was pointed out in reference 1, these influence functions can be found 
either theoretically (refs. 2, 5, and 4) or, if necessary, experimentally. 

If an elastic axis (defined as a line along which loads can be placed 
without producing significant twist anywhere) exists, this equation can 
be simplified to be 
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e(y) 



GM(y;»Tl) 0,(7) d-Tl 


( 2 ) 


The quantity Q(y) is the section torque about the elastic axis and is 
given by 


Q(y) = M(y) - e(y) c(y) L(y) 


( 3 ) 


in -which e(y) is the distance measured forward from the y-axis -to the 
elastic axis, expressed as a fraction of the local chord. 


Aerodynamic Loads 

The section lift and pitching moment can be expressed in coefficient 
form aa 


L(y) = q. c(y) c^(y) 
M(y) = q c2(y) cm(y) 


(i*-) 


By £LSS-uming linearity of the aerodynamics, the loading coefficients for 
steady roll can be writ-ten as 


C;j(y) 

Cm(y) 


Cmci^(y) ^ + CjngCy) 6 


(5) 


The principal way in which this analysis departs from that of reference 1 
is in the manner of obtaining the loads due to struct-ural deformation 
(the first term on the right-hand side of eqs. ( 5 ))* ‘ttie method of 
reference 1 , the loads due to the arbitrary angle-of-attack distributions 
which arise from structural distortion were de-termlned by an exact applica- 
tion of lifting-surface theory; in the present method, these loads are 
calc-ulated approximately as can be seen from the folJLowing development. 
Consider, for example, the meaning of C 7 (y) . This function is actually 

“e 
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the ratio of the section lift coefficient due to structural distortion 
to the angle of twist 9 (y). In general, C7 (y) is dependent, in an 


apparently complex fashion, on the shape of the twist curve; different 
shapes yield different values of this function. Fortunately, however, 
the value of Ct (y) is relatively insensitive to changes in the shape 




of 6 (y)» This fact suggests that 


“e 


can be adeq_uately approximated 


by calculating it for an angle -of -attack distribution that is reasonably 
close to the expected actual mode shape and then by considering this 
qiiantity to be fixed with respect to changes in the mode shape. This 
procedure of using an effective lift-curve slope, which has been used in 
the past by many investigators (for example, see ref. 5)^ obviously a 3 J.ows 
a considerable simplification of aeroelastic analyses. 


For the rolling problem, the section rolling derivatives c^ (y) 
and Cmp(y) J which must be determined for i:ise in equation (5)^ can also 
serve as a convenient basis for determining Ci (y) and Cjjj (y) . The 

Ce Ctg 

angle -of -attack distribution that yields these coefficients , 


<x(y) 


al + y 
( 1 + a) 2 


( 6 ) 


is, for the present aerodynamic purposes, a fair approximation to the 
actual expected mode shape provided that a, the ratio of the fuselage 
radius to the exposed semispan of the wing, is not too large. Thus, the 
following expressions for used: 


ci^(y) = 


( 1 + a) 2 
aZ + y 


•s 

czp(y) 




( 1 + a) 2 
'al + y 


cmp(y) 

/ 


( 7 ) 


It follows from the foregoing development that the only aerodynamic 
information necessary for the aeroelastic analysis of the rolling problem 
is the section rolling and aileron derivatives. At supersonic speeds. 
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■bhese quantities are readily obtainable for most reasonable configurations. 
For rectangular wings, a rather complete derivation of the section lift 
and pitching-moment coefficients due to rolling and aileron deflection 
is included in reference 1 . Lift and pitching -moment distributions due 
to roll can be obtained from references 6 and 7 for a wide variety of plan 
forms; the lift distributions are given directly, and the pitching -moment 
distributions can be found by proper integration of pressure distributions. 
Aileron loads can be found by methods such as those illustrated in ref- 
erence 8; in some cases, two -dimensional strip theory should be adequate. 

For siib sonic speeds, no such con^ilete coverage has been made. In 
the first place, all the theoretical approaches are approximate to some 
extent. In addition, not nearly so large a variety of plan -form shapes 
has been analyzed. However, papers such as reference 9 afford a consider- 
able amount of help in finding the desired aerodynamic derivatives. 


Aeroelastlc Equations 

If the expressions for the loads (eqs. (ij-) and (5)) are Bi±>stituted 
into the equllibriimi equation (eq. (l)), the following aeroelastlc equa- 
tion results ; 


e(y) = <i 


0 L 




e(ri) dTi + 


, I r 


pb 

2 V Jo L 


GL(y^il) (ti) + c2 Gj4(y,Ti) c^(t)) 


dTi + 


qS 



c QL(y;il) GM(y,Ti) 


-ms 


in) 


dT] 


(8) 


For steady roll, the total rolling moment m'ust be zero. Thus, 


0 = / c(aZ + T]) 


9(11) + ^ c^p(Ti) + 5 Cjg(Tl) 


dll (9) 
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SOLUTION OF AEROELASTIC EQUATIONS 
Galerkin-Type Procedure 


If, for a particular configuration, the values of q and 5 are 
given, equations (8) and (9) can be solved simultaneously to yield the 
values of twist 0 and rate of roll pb/2V for the elastic aircraft. 

In reference 1, these equations were solved by a collocation procedure 
that involved the solution of high-order matrix equations. A method 
that is considerably simpler (which takes the form of the Galerkin method) 
is used herein. The solution proceeds as follows: 

Let 


0(y) = k 0q(y) 


( 10 ) 


where 8q(y) is an approximation of the actual expected twist shape. 

Note that, altho\ogh the Galerkin solution generally involves the use of 
a series of such functions, only one term is used for this particular 
application. 

If the approximation for 0(y) (eq. (lO)) is introduced into equa- 
tion (8), and the resulting eq\iation is multiplied by 9q(y) and inte- 
grated over the exposed semispan of the wing, the following equation is 
obtained: 


k r 0l^(y) dy = q,k r f 9i(y) 
Jq Jq ■Jq 


C GL(y;.il) 


0l(n) it dy + 


^ 2V Jo Jq 


l(y)jc 0L(y,il) GM(y>il) Cjmp(Ti) 


dT| dy + 


qb f f 0i(y)rc GL(y»n) OM(y^n) 

J Q J Q L 


( 11 ) 
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Introduction of equation (lO) into equation (9) yields 


0 = k r c(al + q) (q) eqCq) dq + f c(al + q) c^ (q) dq + 
J Q cxg 2V J Q p 


6 / c(aZ + q) (q) dq 


( 12 ) 


0 


Dividing equations (U.) and (12) by 5 yields tvo simultaneous 

equations in two unknowns, k/6 and — /d. The quantity that is actually 

2V/ 

sought is the latterj solving for this quantity gives 


5^/6= - 

2Y 


aCiq - ^(BeCig - BsC^q) 
- q^B0C2p - 


(15) 


where 


A = / 6-, (y) dy 

^0 


(15a) 


Be = / f 0i(y) 


n z 

9i(y) 


C GL(y^T|) QM(y:-n) <qnae(’^^ 


0l(q) dq dy 


1 Z pZ 


Bq= j P3_: 


01 (y) 


c GL(y,q) Cl (q) + GM(y,n) c^Cq) 


C GL(y^’l) °Z 5 (^l) + %(y^q) Cnig(q) 


dq dy 


dq dy 


> (15b) 
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and 


- 

® , 2 - 


2Z^c(l + a) 0 


c(aZ + T^) 02_(ti) dT] 


Cl- = 


^ 22^c(l + a) 0 


c(al + Ti) dTi 


^ (13c) 




1 


22^c(l + a) 0 


c(a2 + ii) cjg(ii) dTi 


The quantities ^Iqj are actmlly the rolling-moment coeffi- 

cients (based on exposed wing area) resulting from a unit mode shape 
a unit pb/2V, and a unit 5, respectively. 


Rolling effectiveness .- The rolling effectiveness ^ is defined as 
the ratio of the rate of roll of the flexible airplane to the rate of 
roll that would occur if the airplane were rigid. The rate of roll for 
the rigid wing is given by 



Therefore^ the rolling effectiveness is 






(li^) 
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Aileron reversal . - The dynamic pressure at aileron reversal can be 
found by setting p ^eq_ual to zero. The result is 


'Irev ~ 


Be 


A 


1 

A 


( 15 ) 


Twist mode shape . - The mode shape should be a reasonably good 

approximation to the actual twist expected. One possible shape that may 
meet this requirement is the twist that would result from the application 
of the aileron loads only. This shape is given by 


1 r 


01 (y) = 


0 L 


C QL(y;.n) %(y,T) C^Cq) dq (l6) 


Simplifications for wings with elastic axis .- When an elastic axis 
exists; the structural equilibrium is expressed by equation (2) rather 
than equation (l). Therefore; the quantities Eg; Bp; and B5 (in 

eqs. (13b) and the definition of the twist mode shape (eq. (16)) can be 
altered by deleting the terms containing Gj^(y;q) in equations (l3b) 

and (16) and by replacing ‘^“^6 with 

Cq^ = Cmg - e Cj^ 


Numerical Evaluation of Integrals 

The Galerkin proced-ure ^just presented involves the calculation of 
a number of integrals. In this section, these Integrals are found 
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numerically by using Simpson's rule -with 10 equal spanwise Intervals. 
Matrix notation is used to facilitate the representation of the integrals. 
The Integration scheme and the matrix notation are similar to those set 
forth in the section entitled "Matrix Operations" in reference 1 . 

In matrix form, the quantities defined by equations (ija), (ijb), 
and (13c) are 
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and 








Cjg - + yj 

2i2c(i + a) 

S 


®i 


C C7 


CZp - [al + yJ 





S 



> 

2 Z^c(l + a) 

1 


p 



( 20 ) 






^^5 2-/ X 

2 l^c(l + a) 

S 

1 . 




j 


In these definitions^ diagonal, and 

column matrices, respectively, made up of the assumed twist shape; 
is an integrating matrix given hy 



( 21 ) 


This particular integrating matrix, of course, has been obtained by 
applying Simpson's rule. Other schemes for numerical integration could 

be used by appropriately modifying [^sj. It is q_uestlonable, however. 
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whether any increase in accuracy obtained by changing the integration 
rule would be worthwhile for this problem. Note that, for simplicity, 
the subscript notation used in reference 1 to denote the location of 
elements in the matrices has been dropped in this paper. 

The mode shape 03_ (eq. (16)) is written in matrix form as 




gl 


c C-, 


% 




Cfflg 


( 22 ) 


Wlien an elastic axis exists, the matrix formulation is considerably 
simplified by the same procedure mentioned previously; that is, in equa- 
tions (19) and (22), the terms involving should be deleted, and 

the column matrices involving coefficients of moments about the y-axis 
shoiild be replaced by coefficients of moment about the elastic axis . 


APPLICATION 

Computational Procedure 


As can be seen from the preceding einalysis, the requisite qimntities 
for determining the aeroelastlc effects on the roll of a particular air- 


craft are the structioral influence coefficients 


and and 


the aerodynamic derivatives c^^, Cnj^, c^^, and c^. (The deriva- 
tives Cl and Cjn^ are given in terms of c^^ and respectively, 

ty eqs. (7).) The influence coefficients are dependent on the structure 
only, but the aerodynamic derivatives vary with Mach number. If a range 
of Mach mmibers is to be covered, therefore, these derivatives must be 
calculated anew for each value of Mq* 


Not only the aerodynamic derivatives but also the mode shape 0q 
(as calculated from eq. (22)) varies with Mq. The variation of the 
derivatives is unavoidable; the variation of Sq can often be circimi' 
vented, however, by calculating 0q for a particular Mach number and 
then by using this same mode shape for the other values of Mq. This 
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process would Involve little loss in accuracy provided that the lift 
and pitching-moment distributions due to aileron deflection do not change 
radically with Mach number. 

With 0p determined, the calculation of the values of the quantities 
appeari 22 g in equations ( 13 ), (l4), and (I 5 ) proceeds in a straight for- 
ward fashion. (The matrix multiplications in eqs. ( 18 ), (l9), and (20) 
should be performed from left to right becatise the variable aerodynamic 
derivatives are the last terms.) With these quantities determined, the 
dynamic pressure at aileron reversal q^ev calculated from equa- 

tion ( 15 ) , and the rolling effectiveness ^ for other values of q can 
then be coa^juted from equation (lil-). This process is repeated for each 
Mach number until the entire range is covered. 


Sample Calculations 

The method derived in this report is applied to two example config- 
urations. Both of these aircraft have two flexible rectangular wings 
mounted diametrically on a long cylindrical fuselage; both aircraft have 
full-span, 0.2-chord, trailing -edge ailerons. Attention is restricted to 
the supersonic-speed regime. 

The two configurations are shown in figure 2 and are designated as 
models 1 and 2; the wings of both models have the same plan-form aspect 
ratio of Z/c = 1.5* The wings differ, however, in that model 1 has a 
rectangular cross section with a thickness ratio of 0.02, whereas model 2 
has an NACA 65 AOO 3 cross-sectional shape. The modeis also differ in 
that the value of a, the ratio of fuselage radius to exposed wing semi- 
span, is 0.2 for model 1 and O .236 for model 2. Both wings were assumed 
to be made of solid aluminum alloy. 

The torsional Influence coefficients for the two models are given 
in table I. These Influence coefficients were obtained from an approxi- 
mate plate theory which is essentially the same as that of reference k; 
however, the analyses of the two models were slightly different: For 

model 1 , the root was assumed to be completely clamped, and no account 
was made for the stiffening flange effect of the bent-up aileron. For 
model 2 , some root flexibility was allowed, and the flange effect of the 
aileron was taken into account approximately. In both cases, the analyses 
indicated the existence of an elastic axis; therefore, only the torsional 

influence coefficients due to torque given in table I. 

As is evident from the preceding discussion, model 1 represents an 
idealized configuration (actually the same as the model considered in 
ref. 1), and model 2 represents a more realistic aircraft. The following 
table summarizes some of the information necessary for analyzing the models 
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Model 1 Model 2 



a 

e 


1.5 

0.2 

0.2 

0 


1.5 

0.2 

0.236 

0.04^ 


No absolute dimensions have been specified because^ as is shovn in ref- 
erence 1, only the ratios are needed to analyze the aeroelastic-rolling 
problem. 


The aerodynamic rolling derivatives were obtained for Mq = 1.108^ 
1.202, 1.338^ 1.667, and 2.848 for which tables are available in ref- 
erence 1. For Illustrative purposes, the values of c^p and Cg^^ for 


the two models at Mq = 1.202 are given in table II for values of 
between 0 and 1.0 in 0.1 increments. The corresponding values of 



and cq^ 


are also Included in this table. 


The aileron derivatives c^^ and c^^ were fomd by assuming that 

two-dimensional theory is adequate for all stations except at the tip 
where the loads are zero. With this assumption, the aileron loads become 


;, == k ^ ^ = 0, 0.1, 0.2, ... 0.9 

''5 p c I 


= 0 


y = 1.0 

I 


(23) 


'15 


- — fl + 2e - — ) = 0, 0.1, 0.2, ... 0.9 

P c V c / I 


( 24 ) 


- 0 


y = 1.0 

I 


where 
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The asstmed twist mode shape was calculated from eq_uation (22) and 
is given in table IH for the two models. The shapes have been normalized 
by dividing by the tip ordinate and apply for n.T 1 Mach numbers as a con- 
sequence of the Invariance in the shape of the assumed aileron derivatives 
given by equation ( 2^1-) . In the calculation of the mode shape from equa- 
tion (22), the simplifications resulting from the existence of an elastic 
axis were employed. These simplifications were also used wherever else 
applicable. 

With these mode shapes, then, the values of A, Bq, Bp, Bg, 

C-jp, and C-jg were computed from equations (l8), (l9)j and (20) for 

each Mach number. These quantities are tabulated in table IV for the 
two models. From these quantities, the dynamic pressure at reversal and 
the rolling effectiveness ^ can be calculated by equations (I 5 ) and (l^), 
respectively. 


EESULTS AND COMPARISONS 


The results of the aileron-reversal calculations for the two models 
are shown by the test-point symbols in figure 5* la this figure the 
results are given in the form of a plot of the pressure ratio at rever- 


sal 


rev 


against Mach nuaiber, where 




is the static pressure at 


reversal. 


Ph = 


2 

7Mo2 


*3-rev 


and Pq is the standard sea- level static pressure, 2,ll6 Ib/sq ft. For 

comparison, the results obtained by the method of reference 1 are also 
shown In the figure. The agreement is seen to be very good, particularly 
at the higher Mach numbers. 


For the two example configurations considered, the aileron-reversal 
results alone provide an adequate test of the accuracy of the method of 
this paper. This fact arises from the virtual linearity of the variation 
of rolling effectiveness 0 with dynamic pressure (or pressure ratio). 
The rolling-effectiveness curves for model 1 are given in reference 1 and 
are almost linear. For model 2, the calculations by the method of ref- 
erence 1 exhibited even better linearity. Similar degrees of linearity 
also resialt from the method contained herein. For the above reasons, no 
rolling-effectiveness plots are included in this paper. 
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CONCLUDING REMARKS 


The simplified method outlined in this paper for the prediction of 
aeroelastic effects on roll is evidently capable of yielding results 
that compare favorably with those of highly refined methods. Although 
the method has been tested for only two conf igiuratlons at supersonic 
speeds , there is no reason to suspect that the agreement for other con- 
figurations at other speeds would be significantly worse. 

The foregoing- discussion is not meant to imply that this method is 
applicable in all cases. For instance^ one of the most worrisome prob- 
lems facing the aeroelastician is that of chordwise distortion of the 
wing^ the effects of chordwise distortion, which often appears in wings 
with very low aspect ratio, are not considered in this paper. In addi- 
tion, the single-mode Galerkin-type approach used herein may not be good 
enough for some configurations; with a highly swept wing having Inboard 
ailerons, for example, the actual twist distribution changes radically 
■with dynamic press-ure, and no single assumed twist mode shape could be 
expected to yield good results over the entire range of dynamic pressure. 
However, such configurations are rarely encountered. Lastly, it is clear 
that any results obtained by this method woiild be only as good as the 
structural and aerodynamic ingredients introduced into the calculations. 
For this reason, the methods of structural and aerodynamic analysis must 
be reliable. In some cases - at transonic speeds, for instance - resort 
would have to be made to experiment to determine parts of the basic 
information. 

Although the attention throughout this report has been confined to 
the rolling problem, the same type of approach could be used for other 
static-aeroelastic problems such as torsional divergence and center-of- 
pressure shift. 

Langley Aeronautical Laboratoiy, 

National Advisory Conanittee for Aeronautics, 

Langley Field, Ta., December 16, 195^. 
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TABUS I 


T0R3I0JTAL IRFLUHfCB COEETICIEaTS FOR 
S5CAHPLB CONFiaURATIONS 

(a) Model Ij — - 5,151 Ib/eq. ft 








o4 







0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0.007838 

0.015492 

0.019693 

0.021999 

0.023266 

0.023564 

0.0A331 

0.024570 

0.024704 

0.024812 


0 

0.013492 

0.042859 

0. 061202 

0.071283 

0.076821 

0.079869 

0.061560 

0.082320 

0.083104 

0.083551 


0 

0.019693 

0. 061202 

0.102652 

0.129058 

0.143435 

O.1313& 

0.133736 

0.138250 

0.139770 

c. 160879 


0 

0.021999 

0.071283 

0.12903a 

0.179996 

0.2WT25 

0.227862 

0.237366 

0.242760 

0.246047 

0.248444 

W-s 

0 

0.023266 

0.076821 

0.143433 

0.210725 

0.266557 

0.300®9 

0.319369 

0.330304 

0.336645 

0.341613 


0 

0.023964 

0.079869 

0.131360 

0.227862 

0.300639 

0.339^35 

0.395763 

0.416260 

0.428743 

0.437852 


0 

0.024351 

0.081560 

0.153736 

0.237366 

0.519369 

0.395765 

0.457123 

0.493419 

0.518749 

0.555766 


0 

0.024570 

0.082520 

0.138230 

0.242760 

0.330304 

0.41626b 

0.495^19 

0.560052 

0.603147 

0.634583 


0 

0.024704 

0.0©104 

0.159T70 

0.246047 

0.336043 

0.428743 

0.518749 

0.603147 

0.676153 

O.7338S1 


0 

0.024802 

0.083531 

0.160879 

0.248444 

0.541615 

0.437852 

0.355766 

0.634383 

0.733861 

0.633335 





(b) Modal 2) 

^ Ib/fl4 ft 

0** 







0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0.049936 

0.071161 

0.082136 

0.086913 

0.08S996 

0.089903 

0.090298 

0.090473 

O.OS0534 

0.090601 


0 

0.071161 

0.1510 88 

0.162377 

0.176000 

0.181953 

0.184510 

0.18)646 

O.1B6144 

0.166373 

0.186508 


0 

0.082136 

0.162377 

0.226297 

0.259266 

0.273642 

0.279504 

0.282638 

0.2®845 

0.284403 

0.284726 


0 

0.086913 

O.176OCD 

0.259256 

0.324466 

0.338133 

0.3726$ 

0.379071 

0.381688 

0.383355 

0.583950 


0 

0.088996 

0.181935 

0.273642 

0.558053 

0.423736 

0.457373 

0.472337 

0.478846 

0.481867 

0.483612 


0 

0.089903 

0. 184518 

0.2T9SO4 

0.572685 

0.457375 

0.923316 

0.557309 

0.372492 

0.379448 

0.563464 


0 

0.090298 

0.185646 

0.282538 

0.379071 

0.472337 

0.557505 

0.625719 

0.658175 

0.674166 

0.683400 


0 

0.090473 

0.185144 

0.2^843 

0.331088 

0.478846 

0.572492 

0.658175 

0.725414 

0.762137 

0.783372 


0 

O.OSO534 

0.186375 

0.284403 

0.383193 

0.480.867 

0.379448 

0.674166 

0.762137 

0.834629 

0.883359 


0 

0.090601 

0.186508 

0.284725 

0.383930 

0.483612 

0.5S3464 _ 

0.685400 

0.783372 

0.883359 

0.9^351 
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TAHLE H 


AERODYNi\MIC CXDEFFICrEMTS FOR Mq = 1.202 



^2p(y/^) 

cqp(yA) 





Model 

1 


0 

-1.000000 

0 

6.000000 

0 

.1 

-1.1^-15569 

-.038835 

5.662275 

.155335 

,2 

-1. 711-9190 

-.105254- 

5.247568 

.515765 

.3 

-2.015077 

-.185793 

4.836184 

.441105 

.1^ 

-2.212072 

-.266137 

4. 424151 

.532274 

.5 

-2.3335^9 

— 544835 

4.000430 

.591149 

.6 

-2.36737^^ 

-.411680 

5.551063 

. 617521 

• 7 

-2.2926I1-1 

-.456080 

3.056844 

.608105 

.8 

-2.070005 

-.461762 

2.484005 

. 554115 

.9 

-1.6oifOl(-5 

-.595541 

1.749870 

.451281 

1.0 

0 

0 

0 

0 



Model 

2 


0 

-1.1A5629 

0.055563 

6.000000 

-0.291000 

.1 

-1.546673 

.036198 

5.689547 

-.133159 

.2 

-1. 865921 

-.014696 

5.289627 

.041662 

.5 

-2. 117657 

-.080833 

4. 883205 

.186398 

.k 

-2.500751 

-.154066 

4.471253 

. 299412 

.5 

-2.4o86o4 

-.227246 

4.044883 

.381624 

.6 

-2.429104 

-.292778 

5.591357 

.452862 

.7 

-2.341355 

-.341157 

3.091764 

.450477 

.8 

-2.105803 

-.358047 

2.512528 

.427169 

.9 

-1.626495 

-.314764 

1.769672 

.342472 

1.0 

0 

0 

0 

0 
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TABLE III 

TWIST MODE SHAPES tBED IN GALERKEN-TYPE SOLUTION 


jh 

9 x for - 

. 

Model 1 

Model 2 

0 

0 

0 

.1 

.057166 

. 167091 

.2 

. 182638 

. 351836 

.3 

.352681 

. 1)-83678 


. 1 <- 8254 o 

.616914 

.5 

.619788 

. 750274 

• 6 

.757187 

. 822255 

.7 

. 832535 

. 895651 

.8 

.905312 

. 944649 

-9 

.959032 

.978455 

1.0 

1.000000 

1.000000 
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TABLE rv 

VALUES OF A, Bg, Bp, Bg, AND 



1.108 0.14-25598 0.087669 -o.o84o86 -0.0914-806 0.522170 -0.526021 0.462000 
1.202 .425398 .056242 -.053058 -.067719 .5^^5l80 -.552164 .530000 
1.338 .425398 .055716 -.055229 -.050789 .510457 -.514456 .247500 
1.667 .423398 .016822 -.015677 -.033859 .420365 -.41744-9 .165000 
2.848 .425398 .004079 -.003882 -.016950 .252679 -.248355 .082500 

Model 2 

1.108 0.510760 0.128009 -0.112045 -0.174174 0.650158 -0.562615 0.472194 
1.202 .510760 .067671 -.058776 -.124410 .649197 -.581582 .537281 
1.338 .510760 .032038 -.028037 -.093507 .603209 -.537452 .252961 
1.667 .510760 .005056 -.005035 -.062205 .488809 -.433484 .168641 
2.848 .510760 -.006155 .005059 -.051102 .288649 -.256515 .084320 
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Pressire 
ratio at 
reversal, 




Mach number, M 


Flgiire 5 .- Variation of pressure ratio at reversal with Mach number. 


NACA-Lingley - J-29-55 - 1000 




